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Abstract 


A lot of molecules have been reported for certain applications. However, 
the need to continuously search for novel materials with more promise is 
very important because of the competitive technological innovations. One 
of the ways to achieve this is to use valid theoretical methods that can 
effectively and accurately predict the properties of existing molecular 
entities and using these methods to predict the properties of hypothetical 
molecules. The compound (4-methylphenyl) (4-methylpiperidin-1-yl) 
methanone (MPMPM) has been reported for its nonlinear optical 
potentials; however, investigating their reactivity indices would lead to the 
understanding of the mechanism behind their suitability as nonlinear 
optical devices. The molecular backbone of MPMPM was altered by 
introducing some substituents that could alter molecular properties. 
MPMPM and its derivatives were, therefore, optimized with density 
functional theory and the time-dependent density functional theory using 
pure and hybrid correlations and a polar basis set, 6-31G(d). The energy 
band gaps of the substituted derivatives were lower than that of MPMPM 
while the dipole moments and hyperpolarizabilities were higher, indicating 
that they could serve as better alternatives for nonlinear optical 
applications. 


) ie) This work is licensed under the Creative Commons Attribution-Non- 
“ame Commercial 4.0 International License. 





Introduction 


Organic compounds have been investigated 
experimentally and theoretically for their use in 
nonlinear optical (NLO) devices due to the 2-bonds 
in their moieties, non-centrosymmetric structure 
and donor-acceptor (D-A) configuration, which all 
facilitate intra-molecular charge transfer [1-5]. 
These materials shield human parts (eyes, skin, etc.) 
and other optical sensors from dangerous laser light 
by attenuating the incident laser beams [1]. Their 
use has become very important because of the 
growing use of intense lasers in the industry, 
military and medical sectors [1, 2]. These materials 
exhibit properties such as_ two/multi-photon 
absorption [3, 4], nonlinear scattering [5, 6], reverse 
saturation absorption [7, 8], etc. One of the ways to 
improve the NLO ability of a known compound is 
to derivatize it, especially by adding an electron 
donor at one end and an electron acceptor at the 
opposite end, this enhances asymmetric charge 
distribution [9-11]. Asymmetric charge distribution 
can be investigated by calculating the tensor 
components via quantum mechanical calculations 
[12-14]. This is used to calculate the dipole 
moments and hyperpolarizabilities, which are used 
to measure their harmonic generation efficiencies. 
For instance, the molecular first hyperpolarizability 
is used as a measure of the second harmonic 
efficiency of a molecule [12]. Density functional 
theory (DFT) methods have been employed while 
probing into molecular properties, it has been used 
in past and in recent times in different areas of 
materials science research [11, 12, 15-21]. In this 
study, 4-Methylphenyl) (4-methylpiperidin-1-yl) 
methanone, (MPMPM;; Fig. 1) was synthesized and 
reported for its NLO ability [22]. This work aimed 
to corroborate the experimental findings with 
theoretical method(s) and use the method(s) to 
predict the NLO properties of derivatized MPMPM 
to investigate the effect of some substituents on the 
NLO responses via DFT and TD-DFT calculations. 
Their NLO responses would be compared to that of 
urea crystals, a standard for NLO materials. 


Materials and Methods 


MPMPM and its derivatives were modeled and 
optimized using DFT with (i) a restricted Kohn- 
Sham _ self-consistent field (SCF) calculation 
(BLYP) and (11) a restricted hybrid Hartree Fock- 
DFT SCF calculation (B3LYP) [23] using Pulay’s 
direct inversion of the iterative subspace (DIS) and 
geometric direct minimization [24] for both 
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correlations and a polar basis set, 6-31G(d) after 
generating their most stable conformers with the 
molecular mechanics force field. These levels of 
theory have been reported to be consistent with 
experimental trends [25]. The energies of the 
frontier molecular orbitals (FMOs), that is, energies 
of the highest occupied molecular orbital, Exomo 
and that of the lowest unoccupied molecular orbital, 
Etumo were calculated and used to get the values of 
the energy band gap, Eg (eq. 1), chemical hardness, 
n (eq. 2), chemical softness, o (eq. 3) and the dipole 
moments (wu) (eq. 4). The molecular first 
hyperpolarizabilities, ( were calculated from 
Kleimann’s reduced 10 components from the 27 
3x3%x3 matrix (eq. 5) [26]. The absorption 
maxima (Aabs), electronic transition energies/optical 
band gaps (Eex) and oscillator strengths (f) were 
calculated using time-dependent density functional 
theory, TD-DFT method. all on Spartan 14 (v112, 
Wavefunction, Inc., Irvine.). 


Es = E_tumo — Exomo eq. | 
_ Eg 

as eq. 2 
ae 

= eq. 3 

w= (uo + ph + phy /2 eq. 4 

2 
p= (Ca + Byy + Be) + (ae + Brey + 
2 2-4 
Byzz) + (Brzz + Buxz + Byyz) 1/2 eq. 5 


O 


CH 
H3C ; 


Fig. 1 4-Methylphenyl) (4-methylpiperidin-1-yl) methanone 
(MPMPM) and its derivatives (different substituents were 
used to replace the methyl on pheny] ring). For example, -CHO 
replaced -CH3 on the phenyl! ring to become 4-CHO derivative 
(Table 1), CN became 4-CN and so on. 


Results and discussion 


Electronic properties 


The electronic properties are presented in Tables 1 
and Table 2 for DFT and TD-DFT methods. All 
substituted derivatives have lower Eg values than 
MPMPM with nitro and cyano acrylic acid 
substituted derivatives having the lowest values, 
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Table 1 Electronic and nonlinear optical properties of MPMPM and its substituted derivatives with B3LYP/6-31G*. 


B (x10-°° esu) 


Molecules Exnomo (eV) Exumo (eV) 
MPMPM -6.19 -0.61 
4-CHO -6.55 -1.99 
4-CH=C(CN)COOH -6.63 -2.76 
4-CN -6.67 -1.79 
4-NO2 -6.73 -2.58 
-6.58 -1.55 


E,(eV) eV) sev") pO) 
5.58 2.19 0.36 3.52 0.27 
4.56 2.28 0.44 3.80 1.08 
3.87 1.94 0.52 4.44 1.45 
4.88 2.44 0.41 5.73 2.64 
4.15 2.08 0.48 5.87 1.44 
5.03 252 0.39 5.06 2.39 


4-NC 


Table 2 Electronic and nonlinear optical properties of MPMPM and its substituted derivatives with BLYP/6-31G%*. 


B (x10-°° esu) 


Molecules Exnomo (eV) Exumo (eV) 
MPMPM -4.82 -1.15 
4-CHO -5.21 -2.50 
4-~-CH=C(CN)COOH -5.31 -3.15 
4-CN -5.32 -2.26 
4-NO2 -5.39 -3.09 
4-NC -5.23 


-2.06 


which could be attributed to the larger molecular 
weight of the substituted derivatives, mesomeric 
effects and/or conjugation length. The substituted 
derivatives are, therefore, more reactive than 
MPMPM. Nitro is a stronger electron-withdrawing 
group than the others (increase mesomeric effect) 
while cyano acrylic acid has the longest z- 
conjugation length and even the largest of all the 
molecules. The n and o values are expected to have 
similar and reverse trends, respectively (eq. 2 and 
eq. 3). Consequently, all substituted derivatives 
have lower n values, with nitro and cyano acrylic 
acid having the lowest, while the highest o values. 
Formation of instantaneous dipoles (tensor 
components) is usually enhanced by improved 
reactivity of molecules (the closer the HOMO- 
LUMO gap) and the uw values of the substituted 
molecules are higher than that of MPMPM (Tables 
1 and 2) and all molecules under investigation. 
MPMPM showed higher u values than urea (1.98 D) 
[27]. 


Hyperpolarizability 


The B value is a measure of the second harmonic 
generation (SHG) efficiency of a NLO molecule. 
These values are higher in the substituted 
derivatives than MPMPM. A similar trend was 
observed in the two correlations used. The B value 
of urea (a standard for organic NLO materials) was 
calculated as 0.6498 x 10°3° esu [12]. From Table 1 
and Table 2, all substituted derivatives have higher 
values than urea, this 1s because of the enhanced 
delocalization of charges brought about by the 
different substituents that reduced the HOMO- 


AE (eV) n (eV) 5 (eV!) p(D) 

3.67 1.84 0.55 3.43 0.25 
2:41 1.36 0.74 4.05 1.14 
2.16 1.08 0.93 4.69 1.43 
3.06 1.53 0.65 5.83 2.65 
2.30 1.15 0.87 399 1.43 
3.17 1.59 0.63 5.24 2.50 

cD) 

= 

© 

> 

Tb 

cD) 

is 

=) 

) 

© 

O 

200 260 320 380 440 


UV spectrum 


Fig. 2 Calculated UV spectra of the MPMPM and its 
substituted derivatives (MPMPM = black; 4-CHO = red; 4- 
CH=C(CN)COOH = green; 4-CN = blue; 4-NO2 = yellow; 4- 
NC = pink) with hybrid correlation. 


Calculated value 





UV spectrum 


Fig. 3 Calculated UV spectra of the MPMPM and its 
substituted derivatives (MPMPM = black; 4-CHO = red; 4- 
CH=C(CN)COOH = green; 4-CN = blue; 4-NO2 = yellow; 4- 
NC = pink) with pure correlation. 
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Fig. 4 Optimized (4-methylphenyl) (4-methylpiperidin-1l-yl) methanone (MPMPM) (a); optimized cyano derivative (b); 
optimized 4-CHO derivative (c); optimized nitro derivative (d); optimized 4-cyano acrylic acid (e); and optimized isocyano 


derivative (f). 


LUMO gap and improved intramolecular charge 
transfer (ICT). The chemical hardness and softness 
values provide information about the molecules’ 
reactivity and electron transportation. All the 
substituted derivatives are softer than MPMPM as a 
result of the lower energy band gap. This 
consequently led to an increase in their dipole 
moment values and hyperpolarizabilities. 


Electronic transitions 


MPMPM was absorbed at 285 nm with BLYP 
correlation, while it was absorbed at 237 nm with 
B3LYP (Fig. 2 & Fig. 3). The absorption of 295 nm 
was ignored because of its zero-oscillator strength. 
The experimental Abs was 296 nm [22]. The band 
gap corresponding to the calculated aps (BLYP) 
was 4.34 eV while the experimental band gap 
obtained was 4.19 eV. The calculated Aabs (BLYP) 


corresponding to 295 nm (4.21 eV) was ignored 
because of zero oscillator strength. The substituted 
derivatives bathochromically shifted the Abs (Fig. 4, 
& Tables S1-S6). Cyano acrylic acid substituted 
derivative shifted the Aabs the most, followed by — 
NOz. Similar results were observed in previous 
studies [28]. The red shift observed in 4- 
CH=C(CN)COOH is attributed to its conjugation 
length and also the strong inductive electron- 
accepting nature of nitro derivative than others [29, 
30]. 

Conclusions 

DFT and TD-DFT calculations were carried out on 
MPMPM and its substituted derivatives to obtain 
the ground state optimized geometry, FMOs 


energies, energy band gap, dipole moments and 
hyperpolarizabilities. The properties calculated 
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were used to explain their NLO properties and 
compared with the values of urea, a standard NLO 
material. The FMOs energies obtained show that the 
substituted derivatives are more reactive than 
MPMPM owing to lower energy band gaps. The 
dipole moments of the molecules were higher than 
that of urea while the hyperpolarizability values of 
the substituted derivatives were higher than that of 
urea in multiple folds. The hyperpolarizability of 
MPMPM and derivatives depends on the electronic 
Structure of the molecule, as seen from the lower 
energy band gap, which increased charge transfer 
from the donor to the acceptor parts of the 
molecules. It is also evident that the 
hyperpolarizability was affected by the energy band 
gap. The substituted derivatives could be used as 
precursors in the development of better NLO 
materials. In this study, only the electronic 
components of the hyperpolarizability values were 
considered. The effect of field strength on nuclear 
positions was not considered. 
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